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Qv^ , Different modern pliase sliift equivalent NN potentials are tested by evaluat- 

ing the partial wave decomposition of the kinetic and potential energy of the 
deuteron. Significant differences are found, which are traced back to the ma- 
Q ' trix elements of the potentials at medium and large momenta. The influence 

2 ■ of the localisation of the one-pion-exchange contribution to these potentials 

is analyzed in detail. 
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During the past few years, considerable progress has been made in constructing reahstic 
models for the nucleon-nucleon (NN) interaction. Several of these models describe an iden- 
tical data base of NN scattering data with a x^ per datum ^ 1 |l|-§[, meaning in turn that 
the on-shell matrix elements of the NN transition matrix T are essentially equal. This does, 
however, not imply that the models for the NN interaction underlying these descriptions are 
identical. Moreover, the off-shell properties of each potential may be rather different. All 
models for the NN interaction V include a one-pion exchange (OPE) term, using essentially 
the same nNN coupling constant, and account for the difference between the masses of the 
charged (7r±) and neutral (ttq) pion. However, even this long range part of the NN interac- 
tion, which is believed to be well understood, is treated quite differently in these models. As 
we will discuss below, the local approximation for the OPE contribution, employed in the 
non-relativistic Nijmegen []I[ and Argonne potentials, yields significantly different contri- 
butions to the NN potential V as compared to the one derived in the relativistic framework 
underlying the Bonn potential |^. The description of the short-range part is also differ- 
ent in these models. The NN potential Nijm-II [Q is a purely local potential in the sense 
that it uses the local form of the OPE potential for the long-range part and parametrizes 
the contributions of medium and short-range in terms of local functions (depending only 
on the relative displacement between the two interacting nucleons) multiplied by a set of 
spin-isospin operators. The same is true for the Argonne V^ig potential 0. The NN poten- 
tial denoted by Nijm-I |jl]] uses also the local form of OPE but includes a p^ term in the 
medium- and short-range central-force (see Eq. (13) of Ref. [|l|) which may be interpreted as 
a non-local contribution to the central force. The CD-Bonn potential is based consistently 
upon relativistic meson field theory |^. Meson-exchange Feynman diagrams are typically 
nonlocal expressions that are represented in momentum-space in analytic form. It has been 
shown 1^ that ignoring the non-localities in the OPE part leads to a larger tensor compo- 
nent in the bare potential. This is the origin of the fact that the CD-Bonn potential yields a 
smaller D-state probability for the deuteron (Pq = 4.83%) as compared to other potentials 
{Pd ~ 5.6%). 

The aim of this note is to explore characteristic features of the various phase-shift equiv- 
alent NN potentials depending on the model constraints pertinent to the various potentials, 
such as the non-locality of the OPE part of the NN interaction as well as in the exchange 
terms of the heavier mesons. We have here chosen the simplest possible system to study 
these differences, namely the deuteron. Investigations along these lines, using the above- 
mentioned NN interactions, have been made by studying the properties of the triton [Q, 
the optical potential for nucleon- nucleus scattering 0, and the symmetry energy in nuclear 
matter [[7|. Significant differences between non-local and local NN potentials have been ob- 
served in the study of the symmetry energy. As a side-remark, we ought to state that our 
aim, and that of the above studies, is different from earlier studies of phase-shift equivalent 
NN potentials, like e.g., Sauer et al. |§|-0. In those studies, NN scattering matrices were 
constructed to be phase-shift equivalent from a mathematical point of view, that is, the 
structure of the operators which determined the model for the potentials were kept fixed, 
while potential parameters were allowed to vary. 

By construction, all realistic NN potentials reproduce the experimental value for the 
energy of the deuteron of -2.224 MeV. However, the various contributions to the total 
deuteron energy originating from kinetic energy and potential energy in the ^Si and ^Di 



partial waves of relative motion, 

E = {^s\T\^s) + (^d|T|^d) + (^51^1^5) + (*d|1^|^d) + 2(^5|^|^d) 

= Ts + Td + Vss + Vdd + VsD, (1) 

exhibit quite different results. This can be seen from the numbers listed in the upper part of 
Table |. In that table, we display the various contributions to the deuteron binding energy 
employing the four potentials introduced above. For a comparison we also list in this table 
results obtained from older NN potentials. 

The kinetic energies are significantly larger for the local potentials Vis and Nijm II 
than for the two interaction models CD-Bonn and Nijm I which contain non-local terms. 
The corresponding differences in the S'-wave functions can be seen in Fig. ^ where the 
momentum distribution multiplied by the square of the momentum k is plotted for the S 
wave component of the wave function 

ens{k) = k'mlik). (2) 

Note that, besides a factor 1/M with M the nucleon mass, this function corresponds to the 
integrand for evaluating the kinetic energy Tg. 

Comparing the contributions to the potential energy, one finds large differences particu- 
larly for the tensor contribution Vsd- The dominant part of this tensor contribution should 
originate from the tensor component of the OPE potential. This can also be seen from the 
last three lines of Table | where the expectation value of the OPE part of the NN interaction 
is calculated with the deuteron wave function derived from the Bonn potential. These OPE 
contibutions to Vss^ Vdd and Vsd have been evaluated using three different approximations 
to the OPE potential. The first one denoted by "vr Bonn" in Table | corresponds to the 
relativistic nonlocal OPE of the CD-Bonn approach, in which the matrix elements in the 
"^Si-^Si channel between plane-wave states of momenta k' and k are given by 

(,-H.,.,|,) , -^^./^ /' dcos. f^^V "'"^'-^ '^-t - "'\ (3) 

^ ' ^^' / Att 27r M^ ^ EkEk' J -1 \ A^ + q^ J q^ + ml ^^ 

where q^ denotes the momentum transfer, 

g2 = (k - k')^ = k^ + k'^ - 2kk' cos e, 

Ek = VWTap, (4) 

and m^ is the mass of the pion (note that in order to shorten the notation we ignore in these 
equations the charge dependence of m^ and the nucleon mass M); A stands for the cut-off 
parameter, which was chosen to be A = 1.7 GeV 0. Results for such matrix elements as a 
function of k, keeping k' = 95 MeV/c fixed, are displayed in Fig. ^. If we introduce now the 
nonrelativistic approximation for 

E,E', -M'^ \e + h'\ (5) 

we obtain (ignoring the form factor (A^ — ?ti^)/(A^ + q^) and the relativistic square-root 
factors) 



\ I ss\ /local 47i27iM^J-i \2{q^ + ml) 2) ^' 

This can be viewed as the local approximation to the OPE since the matrix element depends 
on the momentum transfer g^, only. It can easily be transformed into the configuration- 
space representation, resulting in a Yukawa term plus a 5 function, which originates from 
the Fourier transform of the constant 1/2 in Eq. (^. The various steps leading to this 
result are shown in Fig. 0. It is obvious from this figure that all of the steps leading to 
the local expression (^) are not really justified for momenta k around and above 200 MeV, 
a region of relative momenta which is of importance in the deuteron wave function (see 
Fig. |l|). This is true for the matrix elements Vss as well as Vsd- It is remarkable that the 
regularization of the local OPE by the local Gaussian form factor introduced in the Argonne 
^18 potential leads to matrix elements in the S'D-channel which are close to those derived 
from the relativistic expression of the Bonn potential. This is not the case in the SS channel, 
where the removal of the 5 function term is very significant. This comparison of the various 
approximations to the OPE part of the NN interaction demonstrates that even this long 
range part of the NN interaction is by no means settled. The local approximation and the 
regularization by form factors have a significant effect. 

As a next step we would like to discuss the importance of the OPE contribution as 
compared to the total NN interaction in the partial waves relevant for the deuteron. As 
an example we consider the meson-exchange model of the Bonn potential (figures in the 
left column of Fig. |^) and the local Argonne V^is potential considering matrix elements in 
momentum space between "^Si-^Si (upper part) and ^Si-^Di partial waves (lower part). The 
predominant modification of the vr-exchange in the OBE model is caused by the exchange 
of the p-meson. This leads to a considerable reduction of the tensor component (see Vsd 
in the lower part of Fig. |^) and enhancement of the central contribution (see Vss)- The 
effect of other mesons is very weak in the tensor channel and leads to a global attraction for 
the central part of the NN interaction. The modification of the OPE part in the Argonne 
potential is quite different. The tensor part of the NN interaction is reduced much less 
than in the OBE model. In the central SS channel the short range components of the Vi^ 
potential introduce considerable repulsion but lead to a momentum dependence which, apart 
from a repulsive constant, is similar to that of the Bonn potential. 

Finally, we compare in Fig. ^ the momentum space matrix elements of all four phase- 
shift equivalent charge- dependent NN potentials. Substantial differences are observed in 
the central SS channel as well as in the tensor channel SD. The local approximation of 
the OPE contribution to the tensor term is highly suppressed at large momenta (compare 
Fig. 1^). The differences between the local approximation of the OPE tensor term and the 
nonlocal treatment shown in Fig. ^ is much larger than the final differences between the 
various potentials shown in Fig. ^ In the OBE model the final shape of the tensor potential 
is achieved by means of relativistic corrections, form factors, and p exchange. The local 
potentials describe the reduction of the local OPE mainly in terms of local form factors. 
From the remaining differences we see that the requirement of identical NN scattering phase 
shifts provides only a moderate constraint on this tensor component. As a consequence one 
obtains different D-state probabilities in the deuteron wave function (see Table |). It seems 
that the range of possible D-state probabilities has been reduced to some extent by the 



improvements in NN scattering data analysis (as compared to the range produced by the 
older potentials). 

The momentum dependence of the various NN potentials in the central SS channel is 
rather similar. Here the main differences can be described in terms of a constant in momen- 
tum space or a corresponding 6 function in configuration space. It is worth mentioning that 
the two non-local potentials Nijm 1 and CD-Bonn yield matrix elements, which are closer 
to each other as compared to the other two potentials. The same is true for the Vdd matrix 
elements. This may also be the origin for the differences in the high momentum components 
of the wave functions (cf. Fig. |ip. 

In summary, although the modern NN potentials yield the same binding energy for 
the deuteron, there are significant differences in the contributions to both the kinetic and 
potential energy in the various partial waves. Speaking in general terms, these differences 
can be traced back to off-shell differences between the potentials. It is well known that 
the off-shell behavior of NN potentials cannot be pinned down by on-shell data. On the 
other hand, within realistic physical models, not everything is possible that mathematics 
may allow. Moreover, within a given model, pinning down the on-shell point limits the 
range of variation for the off-shell behavior. Therefore, since the four models considered in 
this investigation reproduce the NN scattering data with the perfect x^/datum ^ 1, and 
assuming that the differences between the models cover about the range of diversity that 
there is to realistic physicial models for the NN interaction, one may conclude that the off- 
shell uncertainties revealed in this study cannot be further reduced by objectively verifiable 
aspects. Further constraints for the off-shell behavior of NN potentials, if any, may come 
from our preference for a particular theory or from the many-body problem. 
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TABLES 



Pot. 


Ts 


Td 


Vss 


Vdd 


VSD 


Pd 




[MeV] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 


[%] 


CD-Bonn 


9.79 


5.69 


-4.77 


1.34 


-14.27 


4.83 


Argon Vi8 


11.31 


8.57 


-3.96 


0.77 


-18.94 


5.78 


Nijm I 


9.66 


7.91 


-1.35 


2.37 


-20.82 


5.66 


Nijm II 


12.11 


8.10 


-5.40 


0.59 


-17.63 


5.64 


Bonn A 


10.05 


4.41 


-7.29 


0.39 


-9.78 


4.38 


BonnB 


10.02 


5.62 


-5.39 


1.01 


-13.49 


4.99 


BonnC 


10.13 


7.44 


-1.20 


3.46 


-22.05 


5.62 


Argon Vii 


10.53 


8.68 


-1.83 


1.99 


-21.59 


6.08 


Reid 


12.61 


9.53 


-0.47 


4.55 


-28.45 


6.47 


Urbana 


11.15 


6.44 


-6.91 


-0.19 


-12.73 


5.20 


vr Bonn 






2.32 


1.86 


-21.48 




vr local 






-1.99 


2.18 


-26.76 




vr Argon 






-1.65 


1.60 


-18.28 





TABLE I. Contributions to the kinetic and potential energy of the deuteron originating from 
the '^S\ and '^Di parts of the wave function as defined in Eq. ||. Results are listed for the 
charge-dependent Bonn potential (CD-Bonn [^), the Argonne Vi% [^, the Nijmegen potentials Nijm 
I and Nijm II |l|, the Bonn potentials A,B,C [Q|, the Argonne Vu [11|, the Reid soft-core potential 



[12| and the Urbana potential [|I^. Also listed are the expectation values of the charge-dependent 
OPE part, using the non-local approach from the CD-Bonn potential ('vr Bonn'), the local ap- 
proximation ('vr local'), and the local approximation with inclusion of the Argonne form factor 
('vr Argon'). In the latter three cases, these expectation values are calculated by employing the 
CD-Bonn deuteron wave function. The last column of this Table shows the calculated Z?-state 
probabilities of the deuteron. 
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FIG. 1. Momentum distribution in the 5- wave of the deuteron derived from various NN inter- 
actions. To enhance the high momentum components the distribution has been multipUed by k"^ 
as outhned in Eq. (^). 
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FIG. 2. Plane-wave matrix elements of the one-pion-exchange potential (OPEP) using various 
approximations. As an example, the matrix elements in momentum space (A;'|y|/c) are shown as 
functions of k for a fixed value of k' = 95 MeV/c. The left part of the figure exhibits matrix 
elements for the partial waves ^5*1 -^5*1, while the right part shows the tensor component in the 
^S\-^D\ channel with k refering to the momentum in the ^D\ partial wave. 
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FIG. 3. Plane-wave matrix elements of the one-pion-exchange potential (OPEP, solid lines) are 
compared to those for the total interaction (solid lines with squares) for the CD-Bonn potential 
(left part) and the Argonne Vis potential (right part) in the ^Si-^Si (upper part) and ^Si-^Di 
partial waves (lower part). For the Bonn potential also the result of vr plus p exchange (dashed 
line) is shown. For the definition of the matrix elements see Fig. |2[ 
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FIG. 4. Plane- wave matrix elements of various realistic NN interactions in momentum space 
for the partial wave channels ^Si-^Si, ^Si-'^Di, and ^Di-'^Di. For the definition of the matrix 
elements see Fig. pi. 
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